The Growth of the Earliest Supermassive Black 
Holes and Their Contribution to Reionization 



Zoltan Haiman, Mark Dijkstra, and Andrei Mesinger 

Department of Astronomy, Columbia University, New York, NY 10027, USA 

Abstract. We discuss currently available observational constraints on the reionization 
history of the intergalactic medium (IGM), and the extent to which accreting black 
holes (BHs) can help explain these observations. We show new evidence, based on the 
combined statistics of Lyman a and /3 absorption in quasar spectra, that the IGM 
contains a significant amount of neutral hydrogen, and is experiencing rapid ionization 
at redshift z ~ 6. However, we argue that quasar BHs, even faint ones that are below 
the detection thresholds of existing optical surveys, are unlikely to drive the evolution of 
the neutral fraction around this epoch, because they would over-produce the present- 
day soft X-ray background. On the other hand, the seeds of the z ~ 6 quasar BHs likely 
appeared at much earlier epochs (z ~ 20), and produced hard ionizing radiation by 
accretion. These early BHs are promising candidates to account for the high redshift 
(z ~ 15) ionization implied by the recent cosmic microwave anisotropy data from 
WMAP. Using a model for the growth of BHs by accretion and mergers in a hierarchical 
cosmology, we suggest that the early growth of quasars must include a super- Eddington 
growth phase, and that, although not yet optically identified, the FIRST radio survey 
may have already detected several thousand > 1O 8 M BHs at z > 6. 

1 Black Holes and Reionization 

The recent discovery of the Gunn-Peterson (GP) troughs in the spectra of z > 6 
quasars in the Sloan Digital Sky Survey (SDSS) [5,12,51], has suggested that 
the end of the reionization process occurs at a redshift near z ~ 6. On the other 
hand, the high electron scattering optical depth, r e = 0.17 ± 0.04, measured 
recently by the Wilkinson Microwave Anisotropy Probe (WMAP) experiment 
[46] suggests that ionizing sources were abundant at a much higher redshift, 
z ~ 15. These data imply that the reionization process is extended and complex, 
and is probably driven by more than one population of ionizing sources (see, 
e.g., [16] for a recent review). 

The exact nature of the ionizing sources remains unknown. Natural candi- 
dates to account for the onset of reionization at z ~ 15 are massive, metal-free 
stars that form in the shallow potential wells of the first collapsed dark matter ha- 
los [52,6,20]. The completion of reionization at z ~ 6 could then be accounted for 
by a normal population of less massive stars that form from the metal-enriched 
gas in more massive dark matter halos present at z ~ 6. 

The most natural alternative cause for reionization is the ionizing radia- 
tion produced by gas accretion onto an early population of black holes ( "mini- 
quasars"; [21]). The ionizing emissivity of the known population of quasars di- 
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minishes rapidly beyond z ^ 3, and bright quasars are unlikely to contribute sig- 
nificantly to the ionizing background at z 5 [44,18]. However, if low luminosity, 
yet undetected miniquasars are present in large numbers, they could dominate 
the total ionizing background at z ~ 6 [21]. Recent work, motivated by the 
WMAP results, has emphasized the potential significant contribution to the ion- 
izing background at the earliest epochs (z ~ 15) from accretion onto the seeds 
of would-be supermassive black holes [32,41]. The soft X-rays emitted by these 
sources can partially ionize the IGM early on [39,48]. 

In this contribution, we address the following issues: (1) What is the fraction 
of neutral hydrogen in the IGM at z ~ 6? (2) Can quasar black holes con- 
tribute to reionization either at z ~ 6 or at z ~ 15? (3) How did BHs grow, 
in a cosmological context, starting from early, stellar-mass seeds at z ~ 20? 
(4) Can we detect massive early (z > 6) black holes directly? Numerical state- 
ments throughout this paper assume a background cosmology with parameters 
Q m = 0.27, Q A = 0.73, Q b = 0.044, and h = 0.71, consistent with the recent 
measurements by WMAP [46]. 

2 What is the Neutral Fraction of Hydrogen at z ~ 6? 

The ionization state of the IGM at redshift 6 ^, z <, 7 has been a subject of 
intense study over the past few years. While the WMAP results imply that the 
IGM is significantly ionized out to z ~ 15, several pieces of evidence suggest 
that it has a high neutral fraction at z ~ 6 — 7. 

One argument against a simple model, in which the IGM is ionized at z ~ 15, 
and stays ionized thereafter, comes from the thermal history of the IGM [27,47] . 
The temperature of the IGM, measured from the Lya forest, is quite high at 
z ~ 4, with various groups finding values around T ~ 20, 000K [56,34,43]. As 
long as the universe is reionized before z = 10 and remains highly ionized there- 
after, the IGM reaches an asymptotic thermal state, determined by a compe- 
tition between photoionization heating and adiabatic cooling (the latter being 
due to the expansion of the universe) . Under reasonable assumptions about the 
ionizing spectrum, the IGM then becomes too cold at z = 4 compared to obser- 
vations [27]. Therefore, there must have been significant (order unity) changes 
in fractions of neutral hydrogen and/or helium at 6 < z < 10, and/or singly 
ionized helium at 4 < z < 10. An important caveat to this argument is the 
possible existence of an additional heating mechanism that could raise the IGM 
temperature at z ~ 4. Galactic outflows could heat the IGM, in principle, but 
observations of close pairs of lines of sight in lens systems suggest that the IGM 
is not turbulent on small-scales, arguing against significant stir-up of the IGM 
by winds [40] . The known quasar population is likely driving the reionization of 
helium, Hell— ►Helll, at z ~ 3 (e.g. [26]). If this process starts sufficiently early, 
i.e. at z I; 4, then Hell photoionization heating could explain the high IGM 
temperature. It would be interesting to extend the search for Hell patches that 
do not correlate with HI absorption to z ~ 4 to test this hypothesis. 
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A second argument for a large neutral fraction comes from the rapid redshift 
evolution of the transmission near the redshiftcd Lya wavelength in the spectra 
of distant quasars. The Sloan Digital Sky Survey (SDSS) has detected large re- 
gions with no observable flux in the spectra of several z ~ 6 quasars [5,12,51]. 
The presence of these Gunn-Peterson (GP) troughs by itself only sets a lower 
limit on the volume weighted hydrogen neutral fraction of xhi ^ 10~ 3 [11]. How- 
ever, this strong limit implies a rapid evolution in the ionizing background, by 
nearly an order of magnitude in amplitude, from z = 5.5 to z ~ 6 [8,11,30] (we 
note that the Lyman (3 region of the spectra, which is needed for this conclusion, 
is dismissed in another recent study [45], which therefore reaches different con- 
clusions). Known ionizing populations (quasars and Lyman break galaxies) do 
not evolve this rapidly; comparisons with numerical simulations of cosmological 
reionization (e.g. [15]) suggests that we are, instead, witnessing the end of the 
reionization epoch, with the IGM becoming close to fully neutral at z ~ 7. At 
this epoch, when discrete HII bubbles percolate, the mean-free-path of ionizing 
photons can evolve very rapidly (e.g. [22]) and could explain the steep evolution 
of the background flux. 

However, perhaps the strongest argument for a large neutral fraction comes 
from the presence of the cosmic Stromgren spheres surrounding high-z quasars. 
If indeed the intergalactic hydrogen is largely neutral at z ~ 6, then quasars 
at this redshift should be surrounded by large ionized (HII) regions, which will 
strongly modify their absorption spectra [31,7]. Recent work has shown that the 
damping wing of absorption by neutral hydrogen outside the HII region imprints 
a feature that is statistically measurable in a sample of <~ 10 bright quasars 
without any additional assumptions [38] . A single quasar spectrum suffices if the 
size of the Stromgren sphere is constrained independently [38,37]. In addition, 
with a modest restriction (lower limit) on the age of the source, the size of the 
HII region itself can be used to place stringent limits on the neutral fraction of 
the ambient IGM [54]. 

To elaborate on these last arguments, based on the quasar's HII region, in 
Figure 1 we illustrate a model for the optical depth to Lyman a absorption as 
a function of wavelength towards a zq = 6.28 quasar, embedded in a neutral 
medium (xhi = 1), but surrounded by a Stromgren sphere with a comoving 
radius of Rs = 44 Mpc. Around bright quasars, such as those recently discovered 
[10,12] at z ~ 6, the proper radius of such Stromgren spheres is expected to be 
R s « 7.7x HI 1/3 (N Q /6.5 x 10 57 s" 1 ) 1 ^ (t Q /2xl0 7 yr) 1 ^ [(1 + z Q )/7.28]- 1 Mpc 
[31,7]. Here xhi is the volume averaged neutral fraction of hydrogen outside 
the Stromgren sphere and Nq, £q, and zq are the quasar's production rate of 
ionizing photons, age, and redshift. The fiducial values are those estimated for the 
z = 6.28 quasar J1030+0524 [19,54]. The mock spectrum shown in Figure 1 was 
created by computing the Lyman a opacity from a hydrodynamical simulation 
(kindly provided by R. Cen; the analysis procedure is described in [37]). The 
optical depth at a given observed wavelength, A Q b s , can be written as the sum of 
contributions from inside (tr) and outside (td) the Stromg ren sphere, TLy a — 
tr + t_d. The residual neutral hydrogen inside the Stromgren sphere at redshift 
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Fig. 1. Model from a hydrodynamical simulation for the optical depth contributions 
from within (tr) and from outside (td) the local ionized region for a typical line of sight 
towards a zq — 6.28 quasar embedded in a fully neutral, smooth IGM , but surrounded 
by a local HII region of (comoving) radius Rs = 44 Mpc. The dashed curve corresponds 
to to, and the solid curve corresponds to tr. The total Lyman a optical depth is the 
sum of these two contributions, TL ya = tr + to- The dashed-dotted lines demarcate 
the three wavelength regions used for our analysis described in the text. For reference, 
the redshifted Lyman a wavelength is at 8852 A, far to the right off the plot. 



z < zq resonantly attenuates the quasar's flux at wavelengths around A Q (1 + z), 
where X a = 1215.67 A is the rest-frame wavelength of the Lyman a line center. 
As a result, tr is a fluctuating function of wavelength (solid curve), reflecting 
the density fluctuations in the surrounding gas. In contrast, the damping wing 
of the absorption, t/j, is a smooth function (dashed curve), because its value is 
averaged over many density fluctuations. As the figure shows, the damping wing 
of the absorption from the neutral universe extends into wavelengths A Q b s ^ 8720 
A, and can add significantly to the total optical depth in this region. 

The sharp rise in td at wavelengths A Q b s ^ 8720 A is a unique feature 
of the boundary of the HII region, and corresponds to absorption of photons 
rcdshifting into resonance outside of the Stromgren sphere. The detection of this 
feature has been regarded as challenging: since the quasar's flux is attenuated 
by a factor of exp(— TL yQ ), an exceedingly large dynamical range is required 
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Fig. 2. The Keck ESI spectrum of the z = 6.28 quasar SDSS J1030+0524 [51], in units 
of 1CP 18 erg s -1 cm" 2 A -1 , and including uncorrelated la errors. The upper and lower 
panels show the regions of the spectrum corresponding to Lyman a and (3 absorption in 
the redshift range 6.15 < z < 6.28. The Lyman j3 cross section is smaller than Lyman 
a, and a significantly higher column of neutral hydrogen is required to block the flux 
in the Lyman (3 region than in the Lyman a region. This spectrum therefore requires 
a steep increase in the opacity over the narrow range from z = 6.20 to z = 6.16, 
explainable only by a GP damping wing, as shown by the dashed curve in Figure 1. 



in the corresponding flux measurements. However, simultaneously considering 
the measured absorption in two or more hydrogen Lyman lines can provide the 
dynamical range required to detect this feature [37] . 

In particular, we modeled broad features of the Lyman a and Lyman regions 
of the absorption spectrum of the z = 6.28 quasar SDSS J1030+0524. The ob- 
servational input to our analysis is the deepest available absorption spectrum 
of SDSS J1030+0524 [51], shown in Figure 2. The spectrum exhibits a strong 
Lyman a Gunn-Peterson (GP) trough, with no detectable flux between wave- 
lengths corresponding to rcdshifts 5.97 < z < 6.20, as well as a somewhat nar- 
rower Lyman trough between 5.97 < z < 6.18 [51], as shown in Figure 2. The 
flux detection threshold in the Lyman a and Lyman regions of this spectrum 
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correspond to Lyman a optical depths of Ti im ( L y a ) ~ 6.3 and T lim ( Lyi g) w 22.8 
respectively 1 

To summarize these constraints, we have divided the spectrum into three 
regions, shown in Figure 1. In Region 1, with A bs > 8752.5 A, the detection 
of flux corresponds to the upper limit on the optical depth T^ ya < 6.3. Region 
2, extending from 8725.8 A < A bs < 8752.5 A, is inside the Lyman a trough, 
but outside the Lyman (3 trough. Throughout this region, the data requires 6.3 
Ss T hya Sa 22.8. Region 3, with A Q b s < 8725.8 A has a lower limit TL yQ > 22.8. 
As defined, each of these three regions contains approximately eight pixels. 

We modeled [37] the absorption spectrum, attempting to match these gross 
observed features. We utilized a hydrodynamical simulation that describes the 
density distribution surrounding the source quasar at z = 6.28. We extracted 
density and velocity information from 100 randomly chosen lines of sight (LOSs) 
through the simulation box. Along each line of sight (LOS), we computed the 
Lyman a absorption as a function of wavelength. The size of the ionized region 
(Rs) and the fraction of neutral hydrogen outside it (xhi), and the quasar's 
ionizing luminosity, L- mn , were free parameters. Note that changing Rs moves 
the dashed (jd) curve in Figure 1 left and right, while changing xhi moves 
it up and down; changing L; on moves the solid (tr) curve up and down. We 
evaluated tr and tjj for each LOS, and for each point in a three-dimensional 
parameter space of i?s, xhi, and L lon , we computed the fraction of the LOSs 
that were acceptable descriptions of the spectrum of SDSS J1030+0524, based 
on the criteria defined above. 

The procedure outlined above turns out to provide tight constrains on all 
three of our free parameters simultaneously. In particular, we find the allowed 
range for the radius of the Stromgren sphere to be 42 Mpc < Rs < 47 Mpc, and 
a ~ 1 a lower limit on the neutral fraction of cchi ^5 0.17. These results can be 
interpreted as follows. As mentioned previously, the presence of flux in Region 
2 (Thya < 22.8) sets an immediate lower limit on Rs- Region 3, however, yields 
an upper limit on Rs, from the requirement that TL yQ > 22.8 in that region. 
This high optical depth cannot be maintained by tr alone, without violating 
the constraint in Region 1 of TL yQ < 6.3. We note that tight constraints on the 
neutral fraction (xm J> 0.1) can be obtained from the size of the Stromgren 
sphere together with an assumed lower limit on its lifetime [54]. Our direct 
determination of the Stromgren sphere size is only slightly larger than the value 
assumed in [54], lending further credibility to this conclusion. 

Our direct constraint on the neutral hydrogen fraction, xui, on the other 
hand, comes from the presence of flux in the Lyman (3 region of the spectrum 
corresponding to Region 2. Because of fluctuations in the density field (and hence 

1 For our purposes, these optical depths can be taken as rough estimates. Their pre- 
cise values are difficult to calculate, with Tii m (L y /3) especially uncertain [45,30,8,11]. 
However, we have verified that our conclusions below remain unchanged when the 
threshold opacities are varied well in excess of these uncertainties. In particular, con- 
sidering ranges as wide as 5.5 < nimfLya) < 7 and 10 < T lim ( Ly/3 ) < 30 would lead to 
constraints similar to those we derive below. 
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in tr), a strong damping wing is needed to raise TL yQ above 6.3 throughout 
Region 2, while still preserving TL yQ < 6.3 in Region 1. This result is derived 
from the observed sharpness of the boundary of the HII region alone, and relics 
only on the gross density fluctuation statistics from the numerical simulation. 
In particular, it does not rely on any assumption about the mechanism for the 
growth of the HII region. 

Using a large sample of quasars (and/or a sample of gamma-ray burst af- 
terglows with near-IR spectra) at z > 6, it will be possible to use the method 
presented here to locate sharp features in the absorption spectrum from inter- 
vening HII regions, not associated with the target source itself. The Universe 
must have gone through a transition epoch when HII regions, driven into the 
IGM by quasars and galaxies, partially percolated and filled a significant fraction 
of the volume. The detection of the associated sharp features in future quasar 
absorption spectra will provide a direct probe of the the 3D topology of ionized 
regions during this crucial transition epoch (in particular, it should enhance any 
constraint available from either the Lyman a or [3 region alone [14]). 

3 Did Accreting Black Holes Contribute to Reionization? 

The two most natural types of UV sources that could have reionized the IGM 
are stars or accreting black holes. Deciding which of these two sources domi- 
nated the ionization has been studied for over 30 years (e.g. [2]). It has become 
increasingly clear over the past decade that the ionizing emissivity of the known 
population of bright quasars diminishes rapidly beyond z i> 3, and they are 
unlikely to contribute significantly to the ionizing background at z ^ 5 [44,18]. 
This, however, leaves open two possibilities. First, if low-luminosity, yet unde- 
tected miniquasars are present in large numbers, they could still dominate the 
total ionizing background at z ~ 6 [21]. Second, the supermassive black holes at 
z ~ 6 must be assembled from lower-mass seeds which accrete and merge during 
the hierarchical growth of structure. The population of accreting seed BHs can 
contribute to the ionization of the IGM at z ~ 20 [32,41]. 

The above two possibilities can both involve faint BHs that are not individu- 
ally detectable. However, a population of acccting BHs at z 6 would be accom- 
panied by the presence of an early X-ray background. Since the IGM is optically 
thick to photons with energies E below E max = 1.8[(1 + z)/15)]°' 5 x](^ keV, the 
soft X-rays with E <, E mayi would be consumed by neutral hydrogen atoms 
and contribute to reionization. However, the background of harder X-rays would 
rcdshift without absorption and would be observed as a present-day soft X- 
ray background (SXB). Under the hypothesis that accreting BHs are the main 
producers of reionizing photons in the high rcdshift universe, it is a relatively 
straightforward exercise to calculate their contribution to the present-day SXB. 

We assumed for simplicity [9] that the accreting BHs form in a sudden burst 
at redshift z = zq . The total number of BHs was expressed by a normalization 
constant n, defined as the ratio of the total number of ionizing photons emitted 
per unit volume produced by the BH population to the number density of hy- 
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Fig. 3. Constraints on the number of ionizing photons per H atom, r), and the power- 
law index for the slope of the ionizing background, a, based on the intensity of the 
present-day SXB. The quasars are assumed to form at zq = 10 and have a power-law 
spectrum, Fe oc E~ a for E > 13.6 eV. The curves bracket the allowed parameter space 
for the mean (solid line) or maximum (dotted) unaccounted flux in the SXB at ~1 keV. 
The curves shift by +0.1 in a if zq — 15 is assumed. 

drogcn atoms. Full ionization of the IGM at z ~ 6, where recombinations are 
significant, likely requires r\ J> 10 [18]. We find that in order to account for the 
electron scattering optical depth r e ~ 0.17 by partially (prc-)ionizing the IGM 
at z ~ 15, a somewhat smaller r\ is sufficient (see below). 

The spectrum of the ionizing background is a crucial ingredient of the mod- 
cling, and depends on the type of accreting BH that is considered. For luminous 
quasars powered by supermassive black holes, we adopted a composite spectrum 
[42], based on the populations of lower redshifts (0 < z < 5) QSOs. The spectra 
of lower-mass miniquasars, with BHs whose masses are in the range Mbh ~ 10 2-4 
Mq , are likely to be harder. For these sources, we followed [32] and adopted a 
two-component template that includes a multi-temperature accretion disk, and 
a simple power-law emission to mimic a combination of Brcmsstrahlung, syn- 
chrotron, and inverse Compton emission by a non-thermal population of elec- 
trons. 

Finally, we took the unresolved soft X-ray band in the energy range 0.5 — 2.0 
keV to be in the range (0.35 — 1.23) x 10~ 12 erg s _1 cm~ 2 dcg~ 2 . This range was 
obtained [9] from a census of resolved X-ray sources, which we subtracted from 
the total SXB. We included the uncertainties in both the measurement of the 
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total SXB (which was dominant) and in the number of resolved point sources. 
We examined a range of X-ray energies, and found the strongest constraints in 
the 0.5 — 2 kcV band. We note that a recent study [4] of faint X-ray sources 
found a significant population of star-forming galaxies among these sources, 
with a steeply increasing fractional abundance (over AGNs) toward low X-ray 
luminosities. If this trend continues to a flux limit that is only modestly below 
the current point-source detection threshold in the deepest Chandra fields, then 
the SXB would be saturated, strengthening our constraints (leaving less room 
for any additional, high-z quasars). 

We find that models in which z > 6 accreting BHs alone fully reionize the 
universe saturate the unresolved X-ray background at the > 2a level. Pre- 
ionization by miniquasars requires fewer ionizing photons, because only a fraction 
of the hydrogen atoms need to be ionized, the hard X-rays can produce multiple 
secondary ionizations, and the clumping factor is expected to be significantly 
smaller than in the UV-ionization case [18,39]. We find that models in which 
X-rays are assumed to partially ionize the IGM up to x e ~ 0.5 at 6 ^ z <J 20 
are still allowed, but could be constrained by improved future determinations of 
the unresolved component of the SXB. 

As emphasized above, the spectral shape of the putative typical high-z ac- 
creting BH is uncertain; the existing templates, motivated by lower-redshift 
sources, can be considered merely as guides. Figure 3 shows which combinations 
of a (the logarithmic slope of the ionizing spectrum) and r\ are allowed by the 
unaccounted flux in the SXB (the solid and dotted curves cover our inferred 
range of the unresolved SXB). This figure shows that for 77 = 10, a power-law 
shallower than a « 1.2 — 1.4 will saturate the unaccounted flux. For comparison, 
a is in the range— 1.5 ^ a <, — 0.5 for z SS 0.3, and —1.2 <, a <, — 0.6 for 
1 <; z <, 6 for optically selected radio quiet quasars [49] . 

Our constraints derive from the total number of ionizing photons that the 
population as a whole needs to produce to either fully or partially reionize the 
universe. Therefore, our conclusions depend mostly on the assumed spectral 
shape and the required number of ionizing photons per hydrogen atom rj. They 
are independent of the details of the population, such as the luminosity func- 
tion and its evolution with redshift. Future improvements in resolving the SXB, 
improving the limits on the unresolved component by a factor of a few, would 
place stringent constraints on the contribution of z ~ 15 accreting BHs to the 
scattering optical depth measured by WMAP. 

4 How Did Black Holes Grow by Accretion and Mergers? 

Since the seeds of early BHs may have played a role in reionization, it is all the 
more interesting to ask how the earliest BH population was formed, and how it 
evolved. The remnants of metal-free population III stars that form in the first 
collapsed dark halos can serve as the initial <~ 100M Q seeds that later accrete 
and merge together, to give rise to the supermassive BHs making up the quasar 
population at < z < 6, and the remnant BHs found at the centers of local 
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galaxies. Several recent studies [29,35,50,28] have addressed various aspects of 
the evolution of the BH populations, using the underlying merger trees of dark 
matter halos. 

Recent work on the generation of gravitational waves during the coalescence 
of a binary black hole [13,36] has suggested that the binary experiences a typical 
gravitational recoil velocity that may be as large as <; 100 km s _1 . These 
velocities exceed the escape velocity v esc from typical dark matter halos at high- 
rcdshift (z <; 6), and can therefore disrupt the early stages of growth of BHs by 
ejecting the earliest seeds from their host galaxies. BHs can then start growing 
effectively only once the typical dark matter potential wells are sufficiently deep 
to retain the recoiling BHs. 

Relatively little time is available for the growth of fewx 10 9 M SMBHs prior 
to z ~ 6, and their seed BHs must be present as early as z ~ 10 [24]. A model 
in which stellar seed BHs appear in small progenitor DM halos is consistent 
with the presence of a ~ 4 x 10 9 M SMBH at z ~ 10, provided that each 
seed BH can grow at least at the Eddington-limited exponential rate, and that 
the progenitor halos can form seed BHs sufficiently early on [24], in halos with 
velocity dispersions of a ~ 30 km/s. 

We quantified the effect of gravitational recoil on the growth of a few x 10 9 M Q 
black hole [17], by assuming that progenitor holes are ejected from DM halos 
with velocity dispersions a < Vkick/2, and do not contribute to the final BH mass. 
Each halo more massive than this threshold was assumed to host a seed BH that 
accretes at the Eddington rate, and the BHs were assumed to coalesce when their 
parents halos merged. We took, as an example, the SMBH powering the most 
distant SDSS quasar, SDSS 1054+1024 at redshift z = 6.43, with an inferred BH 
mass of ~ 4 x 10 9 M . We find that recoil velocities with Vkick 2i 65 km s _1 must 
occur infrequently, or else this SMBH must have had a phase during which it 
gained mass significantly more rapidly than an Eddington-limited exponential 
growth rate (with a radiative efficiency of ~ 10%) would imply. The super- 
Eddington growth phase can be avoided [55] if seed BHs can form and grow by 
mergers in dark halos with a velocity dispersion as small as a ~ 5 km/s, the halos 
have steep density profiles (increasing the escape velocity [33] from the central 
regions by a factor of J> 5 relative to the naive formula v esc = 2cr) , and the BHs 
that are retained during the mergers of their halos can grow uninterrupted at 
the Eddington rate between their birth and z « 6. 

5 Can We Detect Massive z > 6 BHs Directly? 

A natural question to ask is whether massive BHs at z > 6 can be directly 
detected. While the SDSS has detected a handful of exceptionally bright, few 
x 10 9 Mq black holes (the BH mass is inferred assuming these sources shine at 
the Eddington luminosity), they are likely a "tip of the iceberg", corresponding 
to the rare massive tail of the BH mass function. In X-rays bands, the deepest 
Chandra fields have reached the sensitivity to detect nearly ~ 100 times smaller 
holes ( [23], provided they radiate the Eddington luminosity, with a few percent 
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of their emission in the X-ray bands). However, due to the small size of these 
fields, they have revealed only a handful of plausible candidates [1,3]. 

Detections however, seem promising in the large radio survey, FIRST. We 
used a physically motivated semi-analytic model, based on the mass function of 
dark matter halos, to predict the number of radio-loud quasars as a function of 
redshift and luminosity [25] . Simple models, in which the central BH mass scales 
with the velocity dispersion of its host halo as Mbh oc crfj alo have been previously 
found to be consistent with a number of observations, including the optical and 
X-ray quasar luminosity functions [21,53]. We find that similar models, when 
augmented with an empirical prescription for the radio loudness distribution, 
overpredict the number of faint (~ 10/xJy) radio sources by 1-2 orders of mag- 
nitude. This translates into a more stringent constraint on the low-mass end 
of the quasar black hole mass function than is available from the Hubble and 
Chandra Deep Fields. We interpret this discrepancy as evidence that black holes 
with masses ^ 10 7 M Q are either rare or are not as radio-loud as their more 
massive counterparts. Models that exclude BHs with masses below 10 7 M Q are 
in agreement with the deepest existing radio observations, but still produce a 
significant tail of high-redshift objects. In the l-10GHz bands, at the sensitivity 
of ~ 10/zJy, we find surface densities of ~ 100, ~ 10, and ~ 0.3 dcg~ 2 for sources 
located at z > 6, 10, and 15, respectively. The discovery of these sources with 
instruments such as the Allen Telescope Array (ATA), Extended Very Large Ar- 
ray (EVLA), and the Square Kilometer Array (SKA) would open a new window 
for the study of supermassive BHs at high redshift. We also find surface densi- 
ties of ~ 0.1 deg~ 2 at z > 6 for mJy sources that can be used to study 21 cm 
absorption from the epoch of reionization. These models suggest that, although 
not yet optically identified, the FIRST survey may have already detected several 
thousand such > 10 8 M Q BHs at z > 6. 
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